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f — . ' By choosing the future event horizon as the horizon of the fiat 

FLRW universe, we show that the interacting holographic dark en- 
C*") ' ergy model is able to explain the phantom divide line crossing. We 

show that if one takes the particle event horizon as the horizon of the 
universe, besides describing lu = — 1 crossing, (based on astrophysical 
, data) he is able to determine appropriately the ratio of dark matter 

to dark energy density at transition time. In this approach, after the 
first transition from quintessence to phantom, there is another tran- 
sition from phantom to quintessence phase which avoids the big rip 
singularity. 
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1 Introduction 

Based on astrophysical data, it is believed that the expansion of the universe 
is accelerating PQ . To explain the present inflation one may assume that 70% 
of the universe is composed of a form of energy, dubbed as dark energy [2], 
that permeates all of space and has negative pressure. Many candidates for 
dark energy such as cosmological constant, A (vacuum energy density) [3], 
dynamical exotic scalar fields with negative pressure [1] and so on have been 
introduced in the literature. In the cosmological constant model, the dark 
energy density, p\, remains constant throughout the entire history of the 
universe, while matter density decreases during the expansion. So, in this 
model there must be a rapid transition from matter to dark energy dom- 
inated era. This is in contrast with the astrophysical observations which 
show that dark energy and matter densities are of the same order of mag- 
nitude in the present epoch. This is known as the coincidence problem [5], 
which also arises in dark energy models consisting of non-interacting exotic 
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fields. By considering an appropriate interaction between dark matter and 
dark energy components which converts dark energy into matter, one may 
able to cure this problem [6] . In [7] , in some detail, it was discussed how the 
vacuum energy density couples to the matter fields through matter creation 
pressure. The process couples cosmic vacuum (dark) energy to matter to 
produce future-directed increasingly comparable amplitudes in these fields. 

Some present data seem to favor an evolving dark energy with an equa- 
tion of state parameter (EOS), tua, less than —1 (phantom phase) at present 
epoch from u>d > — 1 (quintessence phase) in the near past [5]. So it may 
be interesting to take into account the possibility that the universe exhibits 
phantom like behavior in the present epoch or in the future, oj = — 1 crossing 
is not allowed in minimally coupled dark energy models [9 J , but in multifield 
models, models with non-minimal coupling between scalar field and gravity, 
and in the framework of the scalar-tensor theory of gravity this transition 
is admissible [10 1 . A question which may be arisen is that why tu^ = — 1 
crossing has been occurred in the present epoch. This can be regarded as 
the second cosmological coincidence problem |11] . Crossing the phantom 
divide line (i.e., crossing u> = —1) may also give rise to big rip singularity 
[12] in a finite future time. 

Based on holographic ideas \1'6\ , one can determine the dark energy 
density in terms of horizon radius of the universe. Following |15j . one can 
show that choosing the particle horizon as the infrared cutoff in the holo- 
graphic dark energy model, leads to a decelerating universe. In [H], the 
future event horizon was examined and it was shown that this choice leads 
to the expected acceleration of the universe. 

To study the acceleration of the universe and the cosmological coinci- 
dence without considering energy exchange between dark matter and dark 
energy, a different approach was proposed in |16j . by assuming that the 
energy density is proportional to the inverse of the area of particle hori- 
zon within a closed universe. Depending on the constant of proportionality, 
either the ensuing inflationary period prevents the particle horizon from 
vanishing, or it may lead to a sequence of big rips as the particle horizon 
repeatedly traverses the closed universe. This model does not require the 
decay of matter, and has a natural reference scale, the radius of curvature. 

In this paper we consider the interacting holographic dark energy model, 
i.e. we assume that the amount of dark energy is proportional to the mass of 
a black hole with the same radius as the event horizon of the universe |14j . 
and besides, assume that there is an interaction between dark matter and 
dark energy (the (interacting) holographic dark energy model was discussed 
vastly in the literature[17J ). We show that by appropriately choosing the pa- 
rameters of interaction between dark energy and cold dark matter, one can 
explain oj = — 1 crossing in flat Friedmann- Lemaitre- Robertson- Walker 
(FLRW) universe. In this paper, lo denotes the EOS parameter of the uni- 
verse. In the first part, the future event horizon is taken as the horizon. 
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We show that this choice, although can describe the phantom divide line 
crossing of the universe, but is inconsistent with the thermodynamics sec- 
ond law for the horizon(thermodynamics of the expanding universe has been 
the subject of several studies [IB]) . In addition, after the transition, the 
universe remains in phantom phase and there is the possibility to encounter 
the big rip singularity. 

Subsequently, we show that if one takes the particle event horizon as the 
horizon of the universe, besides describing u = — 1 crossing, by choosing 
suitable interaction parameters (based on astrophysical data), he is able to 
determine appropriately the ratio of dark matter to dark energy density at 
transition time. In this approach, after the first transition from quintessence 
to phantom, there is another transition from phantom to quintessence phase 
which avoids the big rip singularity. 

We use units h = c = G = kb = l throughout the paper. 

2 preliminaries 

We consider the flat FLRW metric and assume that the universe is composed 
of two perfect fluids, cold (pressureless) dark matter and dark energy. We 
consider exchange of energy between these components, so they do not evolve 
independently: 

Pm + 3Hp m = Q 
p d + 3H(l+uj d )p d = -Q. (1) 

H is the Hubble parameter which in terms of the scale factor a{t) can be 
written as H = a(t)/a(t). The over dot indicates the derivative with respect 
to the comoving time. p m and p d are dark matter and dark energy densities 
respectively. uj d is the equation of state parameter of dark energy. Q is the 
interaction term which may be taken as [6] 

Q = (A m p m + X d p d )H. (2) 

A m and are two numerical constants. Whereas dark matter and dark 
energy stress-energy tensors are not conserved, the total stress-energy tensor 
is conserved: 

p + 3H(1 + uj)p = 0. (3) 

p(= Pm + Pd) is the total fluid density of the universe, in terms of which we 
have 

*■ - 

H = -4ur(l+u)p. (4) 
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The EOS of dark energy component can be written as: 

where 0^ denotes the ratio of dark energy density, pd, to the total density 
p: Qd = pd/ P- Note that < < 1, therefore for to < —1, we have always 
u)d < —1, i.e. if quintessence to phantom phase transition occurs, the dark 
energy component must exhibit the phantom like behavior but the inverse 
is not true: The dark energy component may be phantom like: LOd < — 1 
while the universe remains in quintessence phase, —1 < to < —1/3. 

Using Eqs. JT|), ^ and (H]), one can verify that the evolution equation of 
the ratio of densities of dark matter and dark energy, denoted by r = p m /Pd, 
satisfies 

A/ 



Using 



we find 



r(r + l)H 3w + A m + -2 . (6) 



:= ^ = T^j— > (7) 
p 1 + r 



" 3H (1 - fi d ) 3(1 - fid) 3 ■ U 



One can consider pd as the holographic dark energy density [13] 

where c is a positive numerical constant and L is the infrared cutoff. A 
candidate for this cutoff is the future event horizon defined by 

f°° 1 



In the presence of big rip at t = t s , oo in ()10p must be replaced with t s . 
Using 

R f = HR f -l, (11) 
which follows from (jlOp . and -ff-R/ = c/y/Tt^, we arrive at 

1 2^ 1 li d 
3 3c 3ii lid 

To derive (JEJ) we have also used to = -1 - 2H/(3H 2 ). Eqs.flED and (fT2|) 
may be used to determine u; and in terms of £ld'- 

w = n d - Y^d ~ -3-1 ( 13 ) 
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n d = Hn d ( su + 1 + -Vn~ d j . (14) 

Note that, by definition, Q d lies in (0,1). If at a specific point, fid, = 0, 
(|14p implies that higher derivatives of fid must also be zero at that point 
(denoted by the point of infinitely flatness). By considering that is as an 
analytic function, infinitely flatness may only occur at t —* oo. 



3 Future event horizon and u = — 1 crossing 

In the quintessence phase we have: — 1 < uj < —1/3, while in phantom phase 
uj < — 1. At transition time, we have u; = — 1. So if the transition is allowed 
then oj(fld) + 1 has at least one positive root in (0, 1) and u> is a decreasing 
function of time in the neighborhood of this root. In terms of u defined by 
u = y/fld, u> = — 1 becomes 



u 3 +pu 2 + q = 0, (15) 



where 



c c 
P = ^ + Xd ~ 9=2( A ™ -3 )- ( 16 ) 

In order to have quintessence to phantom phase transition, the cubic equa- 
tion (I15p must have at least one root in (0,1). In addition, at transition 
time we must have Cj < 0. 




a, | _ii . ) . (i7) 

leads to 

— ^(^-Od^)- (18) 

at a; = — 1. Therefore a; < 0, at transition time, is satisfied when 

< u < Minimum{c, }, (19) 

3 

or 

2p 

Maximum{c, — — } < u < 1. (20) 

To go further let us discuss the behavior of the future event horizon at 
transition time, denoted by t = to- in the neighborhood of to, we have 
H(t) = H(t ) + h a (t - to)", where h a = (l/a!)if (a) (t ), and a is the order 
of the first nonzero derivative of H at t Q : = d a H/dt a > 0. Note 

that a is an even positive integer and h a is positive [21]. In a universe 
which will remain in phantom phase, the future event horizon is a non 
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increasing function of time: Rf < [22]. Using (llip , we obtain the following 
expansions for Rf-. 

R f (t) = Rf(t ) (l + ^fj(t - t ) Q+1 ) , (21) 

provided that Rf(to) = 0, and 

R f (t) = R f (t ) + (R f (t )H(t Q ) - 1) (t - t ), (22) 

when Rf(t ) ^ 0. In the case (J2U), we have i?/(t) = R f (t )h a (t - t ) a , 
therefore for t > to, Rf(t) > which is conflict with the fact that in phantom 
era the future event horizon is non increasing. Hence at quintessence to 
phantom transition time we must have Rf ^ which, by considering Rt(t> 
to) < leads to Rf(t = to) = H(to)Rf(to) — 1 < 0, therefore using the 
continuity of Rf we conclude that there exists a neighborhood, N of to, for 
which R f (t G N) < . 

It should be noted that the validity of the above results is contingent on 
the fact that the universe will remain in phantom phase for all future times, 
i.e. if there is another transition from phantom to quintessence phase, we 
may have Rf > 0, even in phantom era. If this situation is allowed, (|15p 
must have two roots in (0, 1). Therefore Rolle's theorem implies that there 
exists a t G (0, 1) such that u(t) = 0. At this time we have u = —2p/3 and 

Co = —Qd /(2c\/^(i) < which follows from (|13f) . Therefore uj is a concave 
function and the (possibly) allowed successive transitions are phantom — ► 
quintessence — > phantom. Hence we conclude that even in the presence of 
two transitions, after the quintessence to phantom transition the system 
remains in phantom phase. On the other hand, by considering HRf = c/u, 
Rf > is satisfied when and only when u < c (see (HID ), therefore following 
the above discussion which implies that at transition time Rf < 0, only the 
choice (]20p is allowed. (I20D also implies that there is a lower bound for dark 
energy density at transition time. 

As an example let us assume Maximum{c, —2p/3} = —2p/3. In this 
case, at transition time we must have —2p/3 < u < 1. In order that the 
transition occurs the cubic polynomial Q(u) := u 3 + pu 2 + q, must have 
at least a positive root in (—2p/3, 1). Following Descartes sign rule, for 
p > 0, q < and p < 0, q < 0, Q(u) has a positive root, while for p < 
0, q > it has either two or no positive roots. Consider the sequence 
D = {u 3 + pu 2 + q, 3u 2 + 2pu, 6u + 2p, 6}. We have 

D(l) = {l+p + q,3 + 2p,6 + 2p,Q} 

£>(-f) = {|jr + <Z,0,-2p,6}. (23) 

By Applying Budan-Fourier theorem and Descartes sign rule (and consid- 
ering < — 2p < 3 which results in 1 +p > 4p 3 /27), we conclude that Q at 
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Figure 1: u as a function of O^, for A m = 11/3, = —10/3 and c = 
1/3 (continuous line) and for A m = —3, Xd = —8 and c = 1/6 (points). 

most has one root in (— 2p/3, 1), provided 

4o 3 

^ + 5<0, l+p + </>0. (24) 

As an illustration, the plot of u is depicted for (p = —1/3, q = — 1/2, c = 1/6) 
corresponding to (A m = — 3, A^ = — 8, c = 1/6) and (p = — l,q = 1/9, c = 
1/3) corresponding to (A m = 11/3, A^ = —10/3 and c = 1/3) in fig.([T]). For 
(p = —l/3,g = ~~ 1/2, c = 1/6), q is negative and a; = — 1 has only one 
root which lies in (— 2p/3, 1). In this case the transition is quintessence — > 
phantom. For (p = —l,q = 1/9, c = 1/3), q is positive and lo = —1 has two 
roots in (0, 1). One of these roots lies in (— 2p/3, 1). Note that in this case 
as we have verified previously, the transitions are phantom — ► quintessence 
— > phantom. 

At the end of this section it is worth to note that if we assign an entropy 
to the future event horizon via 

S = TtR}, (25) 

our previous discussions reveal that the second thermodynamics law for the 
horizon is not respected at least in the transition epoch: S(t ~ to) < 0. Note 
that S(t > to) < in phantom phase. The above entropy may be assumed as 
the entropy of the whole system [14^19] . indeed in the holographic approach, 
all information stored within some volume is represented on the boundary 
of that region. 

In another point of view the entropy attributed to the horizon is consid- 
ered as a measure of our ignorance about what is going behind it. In this 
approach the total entropy of the universe, St, is the sum of the entropy of 
perfect fluids inside the horizon, S m , and the horizon entropy, St = S + S m . 



In this way the validity of the generalized second law (GSL) of thermody- 
namics, i.e St > 0, may be investigated [20]. Recently in [21] it has been 
shown that, at least, at transition epoch GSL is not respected. 

In brief in this section we have shown that although by taking the future 
event horizon as the infrared cutoff, the holographic dark energy model can 
describe u = — 1 crossing. Besides, in this model, the universe will remain in 
phantom phase after the transition(s) and the cosmological evolution may 
be ended by the big rip singularity. 



4 Particle horizon and twice cu = — 1 crossing 

In this section we consider the particle horizon defined by 



a(t>) ' 



as the infrared cutoff and study uj = — 1 crossing in accelerated expanding 
universe. In this situation, (|13|) and (|14|) must be replaced with 

« = - 1 + X l~ Xm n d + l4-^ (27) 

and 

n d = Hn d \3u + 1 - - c yftdj , (28) 

respectively. For u> < —1/3, i.e. for accelerated expanding universe, is a 
decreasing function of comoving time. In terms of p and q defined in (|16|) . 
we have 

2 

lo + 1 = — (u 3 -pu 2 -q). (29) 
3c 

Following Descartes sign rule, the equation uj + 1 may have two positive real 
roots only when p > 0, q < 0. If the equation uj + 1 = has two roots in 
(0, 1), following Rolle's theorem we expect that Co{t) = for a t between the 

roots, uj = occurs at u = 2p/3, and at this point uj = il^ /(2c\/Tld) > 0, 
which shows that uj(t) is convex function. At uj = —1, 

Transition from quintessence (phantom) to phantom (quintessence) occurs 
when uj < (>)0, which following (f3"U|) leads to no > (<)2p/3, where no is the 
root of (f29j) . Based on the above discussion we conclude that in order to 
have two transitions, one from quintessence to phantom and the other from 
phantom to quintessence, the equation (I29D must have two roots. One of 
these roots lies in (0, 2p/3) and the other must be located in (2p/3, 1). The 
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minimum of oj occurs at u = 2p/3 in phantom era: at this point Co = and 
to > 0. 

Consider the sequence D(u) = {u 3 — pu 2 — q,3u 2 — 2pu,6u — 2p, 6}, 
we have D(0) = {-g, 0, -2p, 6}, = {1 - p - q, 3 - 2p, 6 - 2p, 6} and 

D(2p/3) = {— 4p 3 /27— q, 0, 2p, 6}. The discriminant of the cubic polynomial, 
u 3 — pw 2 — g, is: —q{27q + 4p 3 ). If we expect that this polynomial has two 
roots in (0, 1), the third root must also be real (p and q are real) and therefore 
the discriminant must be positive. But as we have seen before, q < 0, which 
results in —4p 3 /27 — q < 0. Therefore there is only one sign change in 
D(2p/3). In D(0), for p > q < there are two sign changes. So if we 
assume 1 — p — q > 0, there is no sign change in -D(l) ( note that 2p/3 < 1 ) 
and following Budan-Fourier theorem we conclude that there is one root in 
(0, 2p/3) and the other root is in (2p/3, 1). In brief if p and q satisfy 

An 3 

0<p<1.5, q<0, ^+^>0, l-p-g>0 (31) 

the transitions quintessence — * phantom — ► quintessence are allowed. In this 
case besides that the big-rip is avoided, the thermodynamics second law for 
the horizon is also respected: S > 0, which follows from 

Rp = HR P + 1 > 0. (32) 

At phantom to quintessence transition time, u has a lower bound: 2p/3, and 
at phantom to quintessence transition time we have < u < 2p/3. In the 
accelerating regime, uj < —1/3, and u is a decreasing function of time. 

It is worth noting that there is no sign change in the sequence D{p) = 
{—q, p 2 , 4p, 6}. So the root of (f29|) . corresponding to the ratio of dark energy 
density to total energy density at quintessence to phantom transition, must 
be restricted to u G (2p/3,p). Thus by choosing an appropriate p < 1, one 
can determine r, in agreement with astrophysical data. 

As an illustration, u is depicted in fig.® for p = 0.8, q = —0.01 cor- 
responding to \ m = 2.98, \d = 3.58, c = 1. The transitions occur in 
u = 0.78(0^ = 0.61) and u = 0.12(0^ = 0.014) for quintessence — ► phantom 
and phantom — ► quintessence respectively. 
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Figure 2: u as a function of f^, for A m = 2.98, = 3.58 and c = 1. 
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